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We study the pion spectrum (and in particular taste-symmetry breaking within it) using HYP- 
smeared valence staggered fermions on the coarse and fine MILC lattices (which have asqtad 
staggered sea quarks). We focus on the dependence on lattice spacing and sea-quark mass. We 
also update our results on source dependence. Our main conclusion is that on the MILC fine 
lattices the appropriate power-counting for SU(3) staggered chiral perturbation theory may have 
discretization errors entering at next-to-leading order rather than at leading-order 
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1. Introduction 

In this note we give an update on our results (most recently presented in Ref. ||l|]) for the 
pion spectrum using valence HYP-smeared staggered fermions. These results are needed as input 
into our parallel calculation of Bk (whose status is described in Ref. [Q]), and more generally pro- 
vide information on the systematic errors that are likely to found with improved staggered actions 
(including the "highly improved staggered quark" [HISQ] action [^). 

Smearing the gauge links in the fermion action is a widely used method to reduce discretization 
and perturbative errors. We use the unimproved staggered action with HYP-smeared [Q] links. One 
conclusion of Ref. [jl|] was that this smearing reduced the dominant discretization error — taste- 
symmetry breaking in the pion spectrum — by a factor 2.5 — 3 compared to quarks improved with 
the asqtad action The resulting taste-breaking was thus reduced to the same level as with the 
HISQ action. Thus, for light quarks, where the additional improvements in the HISQ action are 
less important, HYP-smeared quarks represent a simpler alternative to the HISQ action. 

The reduction in taste-breaking was studied in Ref. ^ on the coarse MILC lattices (with 
a f« 0.125fm). For these lattices it was found that, even after HYP-smearing, the size of 0{a^) 
errors was comparable, for the lightest quark masses, to the 0{p^) terms in chiral perturbation 
theory (;tPT). Thus the standard ~ power-counting of staggered chiral perturbation theory 
(S;^PT) must still be used. It was suggested, however, that on the fine MILC lattices (with a ?a 
0.09 fm) taste-breaking would be reduced to the extent that the appropriate power-counting would 
be ^ p'^, thus simplifying the required fitting. Our new results allow us to study this point. 

Another observation of Ref. ^ was that 0{a^p^) effects were significantly smaller than ex- 
pected. With improved statistics we can now quantify this statement. 

Our final new result concerns the dependence of the pion spectrum on the sea-quark mass. 

2. Scaling behavior of pion spectrum 

In Fig. [I), we show m\ as a function of quark mass for both coarse and fine MILC lattices with 
sea quark mass-ratio of m^/m, = 1/5. Parameters for the numerical study are summarized in Table 



gauge action 


1-loop tadpole-improved Symanzik gluon action 


sea quarks 


Nf = 2+\ Asqtad staggered fermions 


valence quarks 


HYP staggered fermions 


parameters 


MILC fine lattices 


MILC coarse lattices 


a 


0.09fm 


0.12fm 


geometry 


28^ X 96 


20^ X 64 


# of confs 


995 


671 


sea quark masses 


ami = 0.0062, 


ami = 0.01, 




ams = 0.031 


am^ = 0.05 


valence quark masses 


0.003, 0.006,..., 0.030 


0.005, 0.01,..., 0.05 



Table 1: Parameters for the numerical study on scaling violations. 
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(a) MILC coarse lattices {a = 0. 12fm) (b) MILC fine lattices (a = 0.09fm) 



Figure 1: vs. 2m^: the data is obtained using HYP-smeared staggered fermions with cubic wall sources 
and Golterman sink operators [jl|]. The legend indicates the tastes. 

lU — we note that the lightest valence quarks have a mass of m^^^^ /\0. We show only results with 
degenerate valence quarks, and display a linear fit. For clarity, we show only the four tastes with 
the smallest errors (those created by pseudoscalar operators which are local in time [jl|]). 

Taste-breaking is substantially reduced on the fine lattices. To make this quantitative, we 
tabulate results for the splitting in the chiral limit 

A(r)= limK(r)-m2(^5)], (2.1) 

with T the taste, in Table ^. Our results are consistent with the S0{4) symmetry expected at 
0{a^) ^ ^, and we quote the result for the taste with the smallest error. The ratio A(r, fine) /A(r, coarse) 
is about 0.3 for all tastes. For HYP-smeared staggered fermions, taste-symmetry breaking is of or- 
der &{a^a^)} Our results are consistent with this expectation, as shown in Fig. ^. 

In light of these results, what is the appropriate power-counting for S;^PT on the fine lattices? 
When considering kaon properties (e.g. Bk), there are loops containing valence dd pions of all 



Taste (r) 


A(r) [(GeV)2] 


A(Fine) /A(Coarse) 


a = 0.125fm 


a = 0.09fm 




0.0278(6) 


0.0087(3) 


0.314(12) 




0.0540(13) 


0.0168(4) 


0.310(11) 




0.0783(17) 


0.0250(6) 


0.319(11) 


I 


0.1005(84) 


0.0300(31) 


0.299(40) 



Table 2: A(r) in physical units on coarse and fine MILC lattices. 



Smearing removes 0((p-as) taste symmetry breaking terms. 
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Figure 2: A(r) vs. a^a^: is obtained at /i = l/a using four-loop running 



Taste (T) 


slope(r) 


c? = 0.125fm 


a = 0.09fm 


^5 


4.258(11) 


4.436(12) 




4.182(14) 


4.391(12) 


^iA 


4.100(20) 


4.352(13) 


^4 


4.031(23) 


4.314(15) 




4.180(33) 


4.414(33) 


^0- 


4.115(42) 


4.373(35) 




4.052(53) 


4.340(37) 


/ 


3.952(73) 


4.312(40) 



Table 3: Slopes, C2, of the linear fits {r\mji{T) Y = c\ + C2{2r\mq) 



tastes. For our lightest valence quark, the ratio of 0{a^) to 0{p^) contributions is characterized by 

(m2(<§^,)-m2(^5))K(^^v) ^ A(^^v)K(^^v) «0.25. (2.2) 

(We use tensor taste as it is the most numerous and its mass lies near the average of the multiplet.) 
Since ;^PT for kaons is characterized by an expansion parameter of {m^^^ / 1 GeV)^ 25%, it may 
be appropriate to use the power-counting 0{a^) ~ 0{p'^), i.e. to treat discretization errors as of 
next-to-leading order (NLO). This will not, however, be appropriate for calculations of fji using 
SU (2) chiral perturbation theory, where the expansion parameter is smaller. 

Further information on the power counting can be obtained by comparing the slopes of the 
linear fits, differences between which are of 0{a^p^). Our results for these slopes for both lattices 
and all tastes are given in Table |3[ We observe that there are significant splittings between slopes for 
different tastes (the differences are more significant than the errors suggest because of correlations), 
but that these splittings are small. For example, the splitting between slopes for tastes i§5 and i§,4 



4 



Scaling behavior and sea quark dependence ofpion spectrum 



Taegil Bae 




are ^ 4% and 2% on coarse and fine lattices, respectively. On the fine lattices, this is consistent the 
power-counting 0{a^) ~ 0{p^), in which case the splittings between slopes are a NNLO effect. 

We also note that SO (4) breaking in the slopes is much smaller than the splittings allowed by 
50(4), although both are effects of the same order, 0{a^p^) [^. This remains to be understood. 

Finally, we note that the difference between the slopes for the two lattice spacings, which is 
roughly 5%, is of the expected size for a discretization effect (i.e. (aAqcD)^)- This is also seen in 
Fig. ^(a). We note that this difference is comparable to the taste splittings on the coarse lattices, 
indicating that the latter are not significantly enhanced over taste-symmetric discretization errors. 



3. Sea quark dependence of pion masses 

In Table we list the MILC coarse-lattice ensembles we use to study the dependence of the 
pion spectrum on sea quark masses. We find no significant dependence within our errors. This is 



Geometry 


am/ (sea) 


a/Mi (sea) 




24^ X 64 


0.005 


0.05 


3.83 


20^ X 64 


0.007 


0.05 


3.78 


20^ X 64 


0.010 


0.05 


4.49 


20^ X 64 


0.020 


0.05 


6.23 


20^ X 64 


0.030 


0.05 


7.56 



Table 4: Parameters for numerical study of sea quark mass dependence on coarse lattices. 

shown for the Goldstone pion in Fig. |3|(b), and for all tastes in Fig. ^ Note that in Fig. ^, different 
colors represent different tastes while different symbols corresponds to different light sea-quark 
masses. 
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Figure 4: (m^ri)'^ vs. 2m^ri(left) and A{T) vs. flm/(sea)(right) 



Taste 


A{T) [(GeV)2] 




12 p 


Cubic wall [CW] 


Cubic U(l) [CUl] 


> 
CD 

(N 


10 - 




0.0278(6) 


0.0274(5) 


8 - 




0.0540(13) 


0.0535(11) 




^4 


0.0783(17) 


0.0779(24) 


1 

o 


6 - 


^5 


0.0253(33) 


0.0240(49) 




4 - 




0.0500(43) 


0.0486(61) 




2 - 




0.0740(57) 


0.0773(101) 






I 


0.1005(84) 


0.1014(134) 




- 




LT-CW LT-CUl NLT-CW NLT-CUl 



Figure 5: A{T) for cubic wall and cubic U(l) sources. 

The righthand plot in Fig. ^ shows A(r) as a function of the light sea quark mass. There is 
a 2 a drop for am/ sea = 0.005, but only improved statistics will allow us to determine if this is 
significant. 

The independence of the pion spectrum on the sea-quark mass is consistent with the results of 
the MILC collaboration with asqtad quarks [^. We note, however, that other quantities, e.g. fj^ 
and /k, do show some dependence on nisea- 



4. Efficacy of different sources 

Finally, we have updated our results on the relative efficacy of our two sources: "cubic wall" 
and "cubic U(l)" sources. (For their definitions see Ref. ||l|].) Both sources project onto specific 
irreps of the cubic group, but we do not know a-priori what the resulting signal/noise ratio will 
be. Here we extend the study of Ref. [|l|] by comparing the sources on the MILC coarse lattice 
ensemble with ami = 0.01 and am^ = 0.05. 

In Fig. ^, we present A(r) for cubic wall (CW) and U(l) sources (CUl). In this plot, LT 
represents "local-in-time" tastes ({i§,i§5,(§4(§;,i§4}) and NLT represents "non-local-in-time" tastes 
({<§4<§5, i§i,/}). Thus, for example, i§(i§5 and 1^4 1§5 form two different irreps of 5W4 but belong to 
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the same irrep of 50(4). S;^PT predicts that the pion spectrum should respect the S0{4) symmetry 
in the chiral limit, both at LO ||[ 0] and at NLO In Fig. |[ we notice that S0{4) symmetry is 
well respected in A(r) since the data for LT tastes and NLT tastes are consistent with each other 
within statistical uncertainty. 

In the case of LT tastes, the results from both sources have comparable errors. However, for 
the NLT tastes, CW sources are clearly preferred, with errors being smaller by about a factor of 2. 
This is a more significant difference than we observed in Ref. [|l|] on quenched lattices. 

5. Conclusions 

We have observed the expected reduction in taste-symmetry breaking by a factor of f» 3 when 
passing from the coarse to the fine MILC lattices. This implies that, for SU (3) staggered chiral 
perturbation theory, using HYP-smeared valence quarks probably allows one to treat discretization 
errors as a NLO effect. Whether this is really the case can only be determined by detailed fitting, 
such as that we are undertaking for Bk- It will also be interesting to do NLO fits of our pion- 
spectrum data itself, but this requires NLO expressions for pions of all tastes, which are in progress 
but not yet completed. 

Our finding that the pion spectrum is almost independent of the sea-quark mass, which is 
consistent with the MILC asqtad results, simplifies our calculation of Bk, since we need not recal- 
culate the entire spectrum for each MILC fine lattice (although we will check this by calculating 
the Goldstone taste). 

6. Acknowledgment 

C. Jung is supported by the US DOE under contract DE-AC02-98CH10886. The research 
of W. Lee is supported in part by the KICOS grant K20711000014-07A0100-01410, by the KRF 
grants (KRF-2007-313-C00146 and KRF-2008-314-C00062), by the BK21 program, and by the 
US DOE SciDAC-2 program. The work of S. Sharpe is supported in part by the US DOE grant no. 
DE-FG02-96ER40956, and by the US DOE SciDAC-2 program. 

References 

[1] T. Bae, et al, Phys. Rev. D77 (2008) 094508, [arXiv:0801.3000]. 
[2] T. Bae, et al, PoS (LATTICE 2008) 275, [arXiv:0809.1220]. 

[3] E. FoUana et al. [HPQCD Collaboration], Phys. Rev. D 75, 054502 (2007) [ai-Xiv:hep-lat/06 10092]. 

[4] A. Hasenfratz and R Knechtli, Phys. Rev. D64 (2002) 034504, [hep-lat/0 103029]. 

[5] K. Orginos et al., Phys. Rev. D60 (1999) 054503, [hep-lat/9903032]; G.P Lepage, Phys. Rev. D59 
(1999) 074502, [hep-lat/9809157]. 

[6] W. Lee and S.R. Shai-pe, Phys. Rev. D60 (1999) 094503, [hep-lat/9905023]. 

[7] C. Aubin and C. Bernard, Phys. Rev. D68 (2003) 034014, [hep-lat/0304014]. 

[8] S.R. Sharpe and R.S. Van de Water, Phys. Rev. D 71 (2005) 1 14505. 

[9] C. Aubin et al., Phys. Rev. D70 (2004) 1 14501, [hep-lat/0407028]. 



7 



